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Abstract 

Recently the Polish- American collaboration OGLE has reported the observation of 
four possible microlensing events by monitoring, over several months, the bright- 
ness of millions of stars in the region of the galactic bulge. If these events are due 
to microlensing, the most accurate way to get information on the mass of the dark 
compact objects, that acted as gravitational lenses, is to use the method of the 
mass moments. Here I apply this method to the analysis of the events detected 
by OGLE. The average mass turns out to be O.28M , suggesting that the lens 
objects are faint disk stars. The same method applied to the five microlensing 
events detected so far by the EROS and MACHO collaborations, which monitor 
stars in the Large Magellanic Cloud, leads to an average value of O.O8M for the 
dark compact halo objects. 
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1. Introduction 



An important problem in astrophysics is the nature of the non-luminous matter in which 
galaxies are embedded. Its presence is inferred from the shape of the measured rotation 
curves. It is well possible that the dark matter in the halo and in the disk of the galaxies 
is made of "brown dwarfs" or Jupiter-like bodies, which are aggregates of the primordial 
elements: hydrogen and helium, and have masses in the range O(10~ 7 ) < M/M Q < O(10 _1 ) 
(De Riijula et al. 1992). Paczyhski suggested a way to detect such objects in the halo and in 
the disk of our own galaxy using the gravitational lens effect (Paczyhski 1986, 1991; Griest 
et al. 1991). 

Recently the French collaboration EROS (Aubourg et al. 1993) and the American- 
Australian collaboration MACHO (Alcock et al. 1993) reported the possible detection of 
altogether five microlensing events, discovered by monitoring over several years millions of 
stars in the Large Magellanic Cloud (LMC), whereas the Polish-American collaboration 
OGLE (Udalski et al. 1993) has found four microlensing events by monitoring stars located 
in the galactic bulge. 

If these observations are true microlensing events, an important question is to determine 
the mass of the dark compact objects, that acted as gravitational lenses. This knowledge is 
also relevant in order to infer the amount of dark matter in our galaxy and in its halo. The 
most appropriate way to compute the average mass is to use the method of mass moments 
developed by De Riijula et al. (1991). 

Here I apply this method to the four events detected by OGLE by looking at the galactic 
bulge and to the ones of EROS and MACHO in the LMC. I also compute for each event 
the value of the most likely mass of the lens object. Of course the very small number of 
observations at our disposal does not yet allow a precise determination of the mass distri- 
bution, so that the present results are preliminary. Nevertheless, it shows how in practice 
it will be possible to get precise information, as soon as a sufficient number of microlensing 
events will be available, on the mass distribution as well as on what fraction they contribute 
to the total dark mass in the halo or in the disk of our galaxy. In the following I will use 
the formulas and the notation derived in De Riijula et al. (1991), to which I refer for more 
details on the derivation; I refer to the paper by Jetzer & Masso (1994) as well, where the 
first three discovered events in the LMC are discussed. 



2. Microlensing rates 



First I compute the microlensing rate T for an experiment monitoring stars in the galactic 
bulge in Baade's window of galactic coordinates (longitude and latitude): I — 1°, b — — 3.°9. 

Let d be the distance of the massive halo or disk object (MHO) to the line of sight between 
the observer and a star in the galactic bulge, t — the instant of closest approach, and vt 
the MHO velocity in the transverse plane. The magnification A as a function of time is 
calculated using simple geometry, and is given by 
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where 

<r = ~—^f- (2) 



d 2 + v%t 2 



E 

The light curve is a universal function determined by the two parameters d/ Re and vt/ Re- 
Re is the Einstein radius, which is 

~ AGMD ... , 2 . , . , 

it E = — x(l - a) = r E fxx(l - x) , (3) 



with M (respectively jj) the MHO mass (in M & solar mass units) and D (xD) the distance 
from the observer to the source (to the MHO); D = 8.5 kpc is the distance to stars in the 
galactic bulge, and r E = 1-25 x 10 9 km. I use here the definition: T = Re/vt (this is slightly 
different from the definition used in De Riijula et al. (1991)). 

The number density of disk stars per unit mass is given by (Bahcall & Soneira 1980) 

dn dn Q ( Dx I sin 6 1 Dx cos b\ dn Q . . 



dM dM 1 \ 300 pc 3.5 kpc J dM 

where the galactic longitude I = 0° has been adopted and ^ = pd |^ f - > with pd{M) = 
O.O5M pc -3 ; H d can be written as follows 

/ xJ (1-11.7 | sin 

^ ) = exp ( 3^ J' (5) 

using also the fact that cos b ~ 1 for the galactic bulge. 

For the contribution of the halo dark matter in the form of massive objects (MHO) located 
in the disk I use the following distribution for the number density per unit mass dn/dM 

dn TT dn a 2 + R 2 GC dn fn , 

= Mx)-TT7 = o , n2 , n o o — — 717' ( 6 ) 



dM v 'dM a 2 + R 2 GC + D 2 x 2 -2DR GC xcosa dM' 

with dn /dM the local dark mass number density. The local dark mass density is p ~ 

8 x 1O~ 3 M pc~ 3 . It is assumed that dn/dM factorizes in functions of M (or //) and x. The 
galactic core radius is a = 5.6 kpc, whereas Rqc = 8.5 kpc is our distance from the centre 
of the galaxy, and a ~ 4° is the angle between the line of sight and the direction of the 
galactic centre. 

In computing Y one must also take into account the fact that both the source and the 
observer are in motion (Griest 1991; Griest et al. 1991). Relevant are only the velocities 
transverse to the line of sight. The transverse velocity of the micro lensing tube at position 
xD is: Vt(x) = (1 — x)v e ± + xv s ±, and its magnitude is 

v t (x) = \J{1- x) 2 | v Q± | 2 +x 2 | v s± | 2 +2x(l - x) | v Q± \ \v s± \ cos 6 , (7) 

where v s ± and v & ± are the source and the solar velocities transverse to the line of sight and 

9 the angle between them. 

For the velocity distribution of the MHOs or the faint disk stars I consider an isothermal 
spherical model, which in the rest frame of the galaxy is given by 

f(v)dh = ^ e-*<** d\ . (8) 
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Since only the transverse velocities are of relevance cylindrical coordinates can be used and 
the integration made over the velocity component parallel to the line of sight. Moreover, 
due to the velocities of the observer and the source, the value of the transverse velocity gets 
shifted by v t (x). The distribution for the transverse velocity is thus 

f(v T )dv T = e-^ T ~* t)2/v « v T dv T , (9) 

where vh is the velocity dispersion for which I adopt the value vh ~ 30 km s _1 (Paczyhski 
1991). The random velocity of the source stars are again described by an isothermal spherical 
distribution, whose transverse velocity distribution is 

g(v s± )dv s± = — g- e- v ^l v l VsL dv s± , (10) 

7TV D 

where v& = 156 km s _1 is the velocity dispersion (Mihalas & Binney 1981; Paczyhski 1991; 
Griest et al. 1991). 

Taking all the above facts into account, F turns out to be ( De Rujula et al. 1991; Griest 
1991) 

(f°° dn§ \ f^ 71 f°° 
Jo ^ ~dji Jo d ® Jo dVs± 9 ^ Vs± ^ 

£ y/x(l-x) Hi{x) e"" 2 | o °° dy y 2 I (2y V ) , (11) 

where y = T](x, 9, v s ±, v Q ±) = and J is the modified Bessel function of order 0. In the 
limit of stationary observer and source star (v & ± = v s ± = 0), rj = and Jo = 1; Hi means 
either Hd or Hh, uth is related to the minimal experimentally detectable magnification 
Ath = A[u = uth}', v&±_ is ~ 220 km s~ x (more precisely it should be multiplied by cos I, 
where I is the galactic longitude but since 1 — 1°, cos I ~ 1). In computing T one should also 
take into account the limited measurable range for the event duration T, which translates 
into a modification of the integration limits. A fact that can be described by introducing 
an efficiency function e(/x) (De Rujula et al. 1991). However for the range of interest here, 
e(fj,) ~ 1, and thus I will neglect it in the present calculations. 

For an experiment monitoring ]V* stars during a total observation time t Q b s the number 
of expected microlensing events is 

N ev = N* t obs T . (12) 
Assuming a delta-function-type distribution for the masses 

dn p 8{[i — ft) 



eq. (|TTD can be integrated. With N± = 10 6 stars and t obs = 1 year one gets 



(13) 



N - = 2 -° 8 ^ 5 x 10-W S pc-3 "™ • < U > 
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for Hi = Hd, and 



N ev = 0.61 y/ft ( - Po ] u TH , (15) 

V 8 x 1O _3 M pc~ 3 / 



for Hi = Hh. The numerical factor in eq. fll4 ) for N ev as a function of b, the galactic 



latitude, varies between 6.8 for 6 = 0° and 1.9 for b = 5°, whereas the factor for N ev of eq. 
(|T5|) remains practically unchanged. 

3. Most likely mass and mass moments 

The probability P that a microlensing event of duration T and maximum amplification 
Amax be produced by a MHO or a faint disk star of mass /i can be derived starting from eq. 
(]TT|) (Jetzer & Masso 1993) and leads to 

P(fi,T) oc — / c$ / dt> s _L fir(u s ±) / (x(l-x)) 2 Hi(x) e~ v 
1 Jo Jo Jo 



exp [ ) lo [2 V — ) . (16) 

One sees that P(fi, T) = P(fi/T 2 ), and that it does not depend on the value of A max . In 
table 1, P is listed for the four events of OGLE for H^ in table 2 the corresponding values 
for the MACHO and EROS events are given, which are computed using the corresponding 
formula for P(fi, T) as discussed in Jetzer & Masso (1994). 

The normalization of P is arbitrarily chosen such that the maximum of P(hmp,T) = 1, 
and ump is the most probable value. The maximum corresponds to \ir\ / 'v 2 H T 2 ~ 213 for 
Hd and ~ 196 for H^. For the LMC the maximum corresponds to [ir 2 E jv 2 H T 2 = 13.0 (with 
vh = 210 km s _1 and r E = 3.17 x 10 9 km). The 50% confidence interval embraces for the 
mass n approximately the range 1/3/xmp up to 3/j.mp- 

A more systematic way to extract information on the masses is to use the method of 
moments as discussed in De Rujula et al. (1991). The moments < fi m > are given by 

< fi m >= ( d/i^/T , (17) 
J dfi 

and < fi m > is related to < r n >= I] e „ ents T n , with r = (vh/te)T, as constructed from the 
observations by 

<r n >= J dN ev r n = Vu TH 7i (m) < /i m > , (18) 
with m = (n + l)/2 and 

V = ANj obs Dr E v H , (19) 

/■27T roc rl rco 

/ d6 / dv sL g(v sX ) / dx Hi(x) (x(l - x)) m e-" / dy y 3 ~ 2m e~ y I (2r,y) . 
Jo Jo Jo Jo 

(20) 

In table 3 the values of jiirn) are listed for m = 0, 0.5 and 1, which are the ones needed 
here. 
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Eq. ( PU[ ) is also useful for computing the average duration < T > (with T 
jj,x{l — x) ) for a microlensing event defined as 



Re 

Vt 



J dN ev ^ J axil — x) 
<T>=- VtV V -. (21) 

N ev 

Assuming a delta-function-type mass distribution (eq. ([13])), one finds 

<T> l = ^^ J ^- (22) 
v H 7i(0-5) 

For disk stars with Hi = H d one gets 

< T > d ~ 52 days , (23) 

and for MHO with H % = H h 

<T> h ~ 42 days . (24) 

The mean local density of MHO (number per cubic parsec) is < fi° >. The average local 
mass density in MHO is < p 1 > solar masses per cubic parsec. The mean MHO mass, which 
one gets from the four events of OGLE, using the faint disk star distribution Hd, is 

<4> = ±±2L H O.28M . (25) 

If instead one uses the H^ distribution, one recovers practically the same value for the 
average, namely: < /j 1 > / < fi° >~ O.29M . The average value clearly suggests that these 
objects are faint disk stars. Notice that 7i(m) depends on the value of the galactic latitude 
b. However the ratios change much less, in particular 7<j(0)/7d(l) varies by at most 10 to 
15% by varying b from 0° to 5°, whereas 7/i(0)/7h(l) remains practically unchanged. The 
same being true for the ratios 7d(l)/7d(0.5) and 7? i (l)/7/ l (0.5). 

Similarly one can also find the mean MHO mass based on the five events detected by 
EROS and MACHO. The corresponding values for 7(771) are 7(0) = 0.280 and 7(1) = 0.0362 
(De Riijula et al. 1991), thus giving 

<fl n > ~ O.O8M . (26) 
< /i° > 

This value is somewhat lower than the one previously computed of 0.14M Q (Jetzer & Masso 
1994), which was based on only three events. The lower value of the MHO mean mass for 
the halo objects indicates that they are brown dwarfs rather than low-mass stars as seems 
to be the case for the events detected by OGLE. 

Although based on few events, these results are of interest. Clearly, once more data will 
be available it will be possible to get more precise information. In particular, it will be 
possible to determine other important quantities such as the statistical error in eq. ( p5|) or 
eq. (|26|) and the fraction / = Mj/p ~ 0.12 pc 3 < ji 1 > of the local dark mass density (the 
latter one given by po) detected in the form of MHOs. 
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Table 1: Values of the most probable mass /imp in units of M Q as obtained by P(/i,T) 
with Hi = H d and T = ^ for the four microlensing events of OGLE. (v H = 30 km s _1 and 
r E = 1.25 x 10 9 km.) 





1 


2 


3 


4 


T (days) 


25 


45 


10.7 


14 




5.2 x 10~ 2 


9.35 x 10~ 2 


2.22 x 10~ 2 


2.91 x 10~ 2 




0.57 


1.85 


0.105 


0.18 



Table 2: Values of the most probable mass [Imp in M Q units for the five microlensing 
events detected in the LMC (Aj = American-Australian collaboration events (i = 1, 2, 3); F 1 
and F 2 French collaboration events). For the LMC: vh = 210 km s -1 and r E = 3.17x 10 9 km. 





A 1 


A 2 


^3 




F 2 


T (days) 


16.9 


9 


14 


27 


30 




9.67 x 10~ 2 


5.15 x 10~ 2 


8.01 x 10" 2 


1.54 x lO" 1 


1.72 x 10" 1 




0.12 


0.03 


0.08 


0.31 


0.38 



Table 3: Values of 7d(m) (for H d ) and 7/ l (m) (for H h ) for m — 0, 0.5 and 1 valid for b = — 3.°9. 



m 


ld{m) 


7/i ( m ) 





252.7 


390.3 


0.5 


14.7 


27.2 


1 


1.6 


2.4 
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